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Abstract 

We study the charmless rare decays B K*K* within the Perturbative QCD 
picture. We calculate not only factorizable and non-factorizable diagrams, but also 
annihilation ones. Our predictions are the following: The longitudinal polarization 
fraction vary from 75% to 99% depending on channels, the branching ratios are of 
order 10"'^ for B°{B°) K*^K*^ and B^ K*^K*^{K*^), much bigger than that 
for B^[B^) K*+K*~{W~^). The direct CP asymmetry in B^ K*^K*^{K*°) and 
B°{B°) K*+K*- is about -15% and -65% if we choose a(v?2) as 95°. There's no 
direct CPV in B^{B^) — > K*^K*^ decays because of the pure b d penguin topology. 
Our predictions will be tested in the future B experiments. 



1 Introduction 



Exclusive B meson decays, especially B VV modes, have aroused more and more 

interest for both theorists and experimenters. Since it offers an attractive opportunity to get 
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a deep insight into the flavor structure of the Standard Model (SM) and the CP violation 
parameters. But things are not so easy due to the Non-perturbative QCD dynamics. Several 
approaches, which include factorization approach (FA) 0121, QCD improved factorizations 
(QCDF ) 0111, Soft-collinear effective theory (SCET) and Perturbative QCD (PQCD) 
IZIIH] approach, have been developed to solve this problem. PQCD is based on kj- factorization 
theorem [3 QHl ^] while others are most based on collinear factorization ^2]. Besides, 
Sudakov factor and threshold resummation E] have been induced in PQCD to regulate 
the End-point singularities, so the arbitrary cutoffs j3] are no longer necessary. 

In the PQCD framework the hard amplitudes for various topologies of diagrams, including 
factorizable, nonfactorizable and annihilation, are all six-quark amplitudes, while in FA and 
QCDF the leading factorizable diagrams involve four-quark amplitudes. This difference leads 
to a different characteristic scale, the former -y/Am^ (~ 1.5GeV) [ZllHj and the latter m^. 
Therefore, we get a larger Wilson coefficients (Cs-e) associated with the QCD penguin in 
PQCD due to the evaluation of the renormalization group, this means penguin diagrams have 
been enhanced dynamically. Current penguin dominated modes data, such as B ^ Kir J3] 
and B ^ (pK seem to fit well with the PQCD predictions [71IH1ITH]. 

The recent B ipK* data ^3 reveal a large transverse polarization fraction, which 
differs from most theoretical predictions and is considered as a puzzle. This indicates B 
VV modes must be more complicated than we think and needs to be investigated more 
thoroughly. Motivated by this, we study another B —>■ VV mode within the Standard 
Model (SM). B —>■ K*K* decays, which governed hy B K* form factors too, may help 
us to know more about the polarization puzzle, as well as the CKM phase angle a JHl and 
new physics. In the following sections, we will perform B —>■ K*K* decays, which have the 
same topologies with B — > KK [201, within the PQCD framework. Our goal is to find out 
the branching ratios, CP asymmetries, as well as the polarization fractions. 

2 Framework and power counting 

In the B — * K*K* modes, the B meson is heavy and sitting at rest. It decays to two 
light vector mesons with large momenta. Therefore the K* mesons are moving very fast in 
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the rest frame of B meson. In this case, the short distance hard process dominates the decay 
amphtude and Final State Interaction (FSI) may not be important in most of the cases, this 
makes the perturbative QCD applicable. For B^{B^) K*^K*~ decay, because of its small 
branching ratio, FSI may occur through intermediate states pp or K*^K*^, etc. If future 
experiments deviate theoretical prediction largely, it might be an indication of strong FSI 
effects. Here we give only the perturbative picture for experiments to test. 

In PQCD approach, the decay amplitude is factorized into the convolution of the mesons' 
light-cone wave functions, the hard scattering kernel and the Wilson coefficients, which 
stands for the soft, hard and harder dynamics respectively. The transverse momentum was 
introduced so that the endpoint singularity which will break the collinear factorization is 
regulated and the large double logarithm term appears after the integration on the transverse 
momentum, which is then resummed into the Sudukov form factor. The formalism can be 
written as: 



~ j dxidx2dx3bidbib2db2hdbsTr[C{t)^Bixi,bi)^K*ix2,b2)^K*ixs,b3) 

H{x,,h,t)Stix,)e~'% (1) 

where the bi is the conjugate space coordinate of the transverse momentum, which represents 
the transverse interval of the meson, t is the largest energy scale in hard function H, while 
the jet function St{xi) comes from the summation of the double logarithms In^Xj, called 
threshold resummation [3 El, which becomes large near the endpoint. 

We use the effective Hamiltonian for the process B K*K* given by |2I] 

Heff = ^ |k [Ci{p)0,{p) + C2{p)02{p)] - VtY,C,{p)Of\p)^ , (2) 

where the CKM matrix elements Vu = V*^Vub, Vt = V*^Vtb, Ci{p) being the Wilson coeffi- 
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cients, and the operators 



O3 



O 



{diUj)v^A{Ujbi)v-A, O2 = {diUi)v-A{Ujbj)v~A, 

idibi)v-A'^{qjqj)v-A, Oi = idibj)v-A'^iqjqi)v-A, 
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{dibi)v^A'^iqjqj)v+A, = {dibj)v^A'^iqjqi)v-A, 



O7 



-{dibi)v-A'^eg{qjqj)v+A, Os = -{dibj)v^A'^eg{qjqi)v+A, 



O9 



-{dibi)v~-A'^eg{qjqj)v-A, Ow = -{dibj)v-A^eq{qjqi)v-A- 



(3) 




i and j stand for SU (3) color indices. 

Now let's analyze these decay channels topologically. First, it is categorized emission 
and annihilation diagrams. Second, each category can be extracted to 4 diagrams, two 
factorizable and two nonfactorizable, in the leading order. Let's take Fig.l{a) for instance, 
the spectator quark can be attached to each of the quark coming from the 4-quark operators 
with a hard gluon. 

For the B^{B^) ]^*o^*o (^gQ^ys (Fig.l), only the operators O3-10 contribute via 
penguin topology with light quark q = s {diagram a) and via annihilation topology with the 
light quark q = d {diagram b and c) or s {diagram d). It is a pure penguin mode with only 
one kind of CKM element, as a result, it will not generate any difference between and 
B^ decay and hence no direct CP violation. Using the PQCD power counting rules |0, We 
can first predict that the main contribution came from the factorizable parts of the emission 
diagram FLf.^{F stands for factorizable, L stands for longitudinal, e stands for emission and 
4 stands for the operator involved) with a large Wilson coefficient C4 + C3/3 — Cio/2 — C9/6. 
The operator Oq disappear here because the vector meson K*^ can not be produced through 
a {S — P){S + P) operator. But in 5 — > K^K^ decay [201; there isn't such constraint and 
the predicted branching ratio is about three times bigger than ours. Second, the transverse 
parts of the emission diagram(F/ve4 and FTei) are down by a factor rk*{rK* — rnK*/mB) or 
rj^*, then the longitudinal parts(FLe4) dominate this process and give a large longitudinal 
polarization fraction. Third, nonfactorizable amplitudes M, including both emission and 
annihilation diagrams, are suppressed by a power of A/ Mb when compared with factorizable 
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ones. At last, we can forecast the factorizable parts of the annihilation diagrams (c, d) 
counteract separately in most of the cases, to be exactly, -^^03(5) and F^aib) vanish and FTa3{5) 
survive but suppressed by r^*, this makes the emission diagram relatively more important. 
The factorizable parts for the space-like annihilation diagram (6) with operator {S—P){S+P) 
and Wilson coefficient Cg + C5/3 — 6*8/2 — C7/6 do not counteract in any case but is still not 
big enough to play the most important role. It is about 10 times smaller than the emission 
ones after calculation. 




Figure 1: Diagrams for B° K*°K*° 

Things are different for the i?^ —>■ K*^K*^{K*^) dec8iys{Fig.2). The operators Os-io 
contribute via penguin topology with the light quark q = s {diagram a) and via the an- 
nihilation topology with q = u [diagram b), while tree operator Oi^2 also contribute via 
annihilation topology [diagram c). We can see that there are two kinds of CKM elements, 
Vu from tree and Vt from penguin, that will induce weak phase and CP violation. We can 
get diagram 2. (a) and 2.(b) by replacing the d quark in Fig.l (a) and (b) to u quark. It 
makes no difference for power counting and the conclusion we have given doesn't change. 
Diagram (c) is a tree diagram, we have a much bigger Wilson coefficient C2 + Ci/3 for 
the factorizable parts and Ci for the nonfactorizable parts, but at the same time, it is an 
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annihilation diagram, as we have stated, and F^r vanish and Ft is suppressed by r\-t 
in this kind of diagram. After taken account of all these two aspects, we can foresee that 
this diagram will be big but not big enough to increase the branching ratio largely, we also 
believe the transverse parts will play a more important role than that of the former channel. 
Our calculation is consistent with our predictions and will be shown in next section. 




(a) (fe) ^ (c) 

Figure 2: Diagrams for B+ K*+K*° 

We put the diagrams of B^{B^) —>■ K*^K*^ decays in Fig.3. From the topology we know 
that Oa-io contribute via annihilation topology with the light quark q = u {diagram h) 
or s {diagram c) and tree operator Oi^2 contribute via annihilation topology {diagram a). 
CPV occurs in this channels for the same reason we have given for — > K*^K*^. But 
when referring to the branching ratio, it is far different from the former two, cause it's a 
pure annihilation mode and the emission diagram which gives the main contribution to the 
branching ratio of the former two channels no longer exist in this process, so we can expect 
a smaller branching ratio for this channel. Besides, the tree diagram involve Ci + C2/3 for 
factorizable parts and C2 for nonfactorizble parts. As is well-known, Ci + C2/3 is small 
(about 0.1 when t = A.SGeV) but C2 is as big as 1.1 {t = A.SGeV), so the factorizable 
parts can not be so important as the former two. Indeed, we found the nonfactorizable tree 
diagram is the biggest one though it is nonfactorizable suppressed after calculation. All our 
calculations fit well with the predictions and they are shown in section 3. 

Now we are going to extract these decay channels within the PQCD framework. For 
convenience. We adopt the light-cone coordinate system then the four-momentum of 



6 




Figure 3: Diagrams for K*+K*~ 

the B meson and the two K* mesons in the final state can be written as: 

Pi = ^(1,1,0t), 

P2 = ^(l-r|„r|„0T), 

^3 = ^(4„1-4.,0t), (4) 



in which tk* is defined by r|^* = i(l - a/1 - 4M|,/M|) ~ Mj^,/M% < 1. To extract the 
hehcity amphtudes, we parameterize the polarization vectors. The longitudinal polarization 
vector must satisfy the orthogonality and normalization : €21 ■ -P2 = 0, e^L ■ -P3 = 0, and 
e2L^ = esL^ = — 1- Then we can give the manifest form as follows: 



e2L = (1 -r^,,-r^,,OT) 



e3L = ^^(-4.,l-r|.,0T). (5) 
As to the transverse polarization vectors, we can choose the simple form: 

e2T = ^(0,0,lT), 

e3T = ^(0,0,lT). (6) 



The decay width for these channels is : 

(7) 



=L,T 



where Pc is the 3-momentum of the final state meson, and |Pc| = ^(1 — 2r|^,). A^"' is 
the decay amplitude which is decided by QCD dynamics, will be calculated later in PQCD 



approach. The subscript a denotes the hehcity states of the two vector mesons with L(T) 
standing for the longitudinal (transverse) components. After analyzing the Lorentz structure, 
the amplitude can be decomposed into: 



W = MIMl + MlMNtlia = T) ■ e*(a = T) + iMre^upa^^e'',* P,'P^ . (8) 
We can define the longitudinal Hq, transverse H± helicity amphtudes 

Ho = MIMl, = MIMn T MI,V?^Mt, (9) 
where r' = %t^- After the helicity summation, we can deduce that they satisfy the relation 

K" 

= \Ho\^ + \H+\^ + \H_\^ (10) 

(7=L,R 

There is another equivalent set of definition of helicity amplitudes 



Aj_ = ^mI,V?^^Mt, (11) 

with ^ the normalization factor to satisfy 

|AoP+ + = 1, (12) 

where the notations Aq, A\\, A± denote the longitudinal, parallel, and perpendicular polar- 
ization amplitude. 

What is followed is to calculate the matrix elements M.N and Air of the operators 
in the weak Hamiltonian with PQCD approach. We have to admit the light cone wave 
functions of mesons are not calculable in principal in PQCD, but they are universal for all 
the decay channels. So that they can be constraint from the measured other decay channels, 
like B — > Kn and B ^ titt decays etc [3 ISj- For the heavy B meson, we have 

^ '^l + MBh5<pBix,b). (13) 



For longitudinal polarized meson, 
1 



[Mk* ^2L(t>K'{x)+ MkA^) + M^/0^.(X)], (14) 
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and for transverse polarized K* meson, 



In the following concepts, we omit the subscript of the K* meson for simplicity. 

The hard amplitudes are channel dependent, but they are perturbative calculable. The 
amphtudes for 5° k*^k*o ^nd K*°K*'^are written as 

\ i=3 i=3,5 

+ ( E MHe^ + E + E ^Hai ) , (16) 

\i=3,5 i=3 1=4,6 

\ j=3 1=3,5 

+V, Mne. + E-^SL + E ^Ha) ■ (17) 

\j=3,5 i=3 j=4,6 / 

respectively, where the subscript H — L, N,T denotes different helicity amplitudes, and e(a) 
denotes the emission(annihilation) topology. The hard parts for the factorizable amplitudes 
F and for the nonfactorizable amplitudes A4 are derived by contracting the wave function 

to the lowest-order onc-ghion-exchange diagrams. 

The helicity amplitudes and A4~ corresponding to — > K*~^K*^ and B~ — > 
K*~K*^ are written as 

M^H = V:Th-V:Ph (18) 
Mh = V^T„-VtPH (19) 

with 

TH-fsF^^l + Mtl 

and 

Ph = fK^Fue. + fBYl + E -^^e. + E ^HaV 

j=4,6 i=3,5 j=3,5 
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The helicity amplitudes for 5° K*^K* and 5° K*^K* are written as 



Kyir _ Y*rnl _ Y* pi 



v: {ibf^hL + mP^,) - v; y: (fB E + E -^Sl) 

q=u,s \ j=3,5 j=4,6 / 



(20) 



= K + Mf^,) -V^Y: (fB E 41 + E -^2L) , (21) 

q=u,s \ 1=3,5 i=4fi / 

where the F^^- and F^^^ vanish and f!^^- takes from of Eqs. ()58l59p . The detailed formulas 
with polarization A4l, -Mn, and A4t for each diagram are given in the appendix. 



3 Numerical analysis 

For the B meson wave function used in eq. lfT^ . we employ the model [HIHl 

4>b{x) = Nbx'^{1 — xY exp 



1 ^ xMs y 

2 ub 



I )2 7,2 
LObO 



(22) 



where the shape parameter ujb = OAGeV has been constrained in other decay modes. The 
normalization constant Nb = 91.784GeV is related to the B decay constant J'b = 0.19GeV^. 
It is one of the two leading twist B meson wave functions; the other one is power suppressed, 
so we omit its contribution in the leading power analysis [21]. The K* meson distribution 
amplitude up to twist-3 are given by pO'^ with QCD sum rules. 



>K*(X 



X] 



(X 



fT 

Ik* 



■.x{l -x)[l + 0.57(1 - 2x) + 0.07C|/^(1 - 2x 



(23) 



|o.3(l - 2x)[3(l - 2xY + 10(1 - 2x) - 1] + 1.68C]/^(1 - 2x) 



2^J2N, 

+0.06(1 - 2a;)^[5(l - 2xf - 3] + 0.36 {1 - 2(1 - 2s)[l + ln(l - x)]}] ,(24) 



^* {3(1 - 2x) [I + 0.2(1 - 2x) + 0.6(10x2 - 10s + 1 
x) + 0.36[l-6x-21n(l-x)]}. 



2v^ 
-0.12x(l 



(25) 
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= ^^x{l -x)[l + 0.6(1 - 2x) + 0.04C2'/'(1 - 2x)], (26) 
V 2iVc 

fT (o 



= l^tl + (1 - 2^) + 0-44(1 - 2x 

+0.4^2^/^(1 - 2x) + 0.88C]/^(1 - 2x) + 0.48[2x + ln(l - x)]} , (27) 

0^*(x) = ■^^^ {3(1 - 2x)[l + 0.19(1 - 2x) + 0.81(10x^ - lOx + 1)] 

-1.14a;(l -x)+ 0.48[1 - 6a; - 2 ln(l - x)]} , (28) 

where the Gegenbauer polynomials are 

4(0 = ^(3^^-1), (29) 
4(0 = ^(35^^-30^2 + 3), (30) 
4(0 = ^(5^^-1). (31) 

In paper [23], Li has suggested to reanalyze the K* meson distribution amplitude in 
order to solve the polarization puzzle of B (pK*. In that channel, Babar ^H] and Belle 
[T7] have reported a longitudinal polarization fraction(i?i) small to 50%, it is different from 
most theoretical predictions and is considered as a puzzle. Many discussions have been 
given |26| |2Zl I2H1 1211 and among which Hsiang-nan Li argued a smaller B K* form 
factor(74o ~ 0.3), which doesn't contradict any existing data, and hence a new distribution 
amplitude for K* meson. Any how, this assumption need to be justified by experiment and we 
will take the traditional wave function in this letter. If future experiment confirms a smaller 
B K* form factor and those argues, we just replace the wave function and get a smaller 
(about 65%) and smaller branching ratios (about 3 x 10^^) immediately. On the other hand, 
if future experiments find a small Rl and branching ratios for B^{B^) K*^{K*~^)K*^, it 
may be a support for a smaller B —* K* form factor and the validity of PQCD. 

We employ the constants as follows J3]: the Fermi coupling constant Gp = 1.16639 x 
lO"^^^^"^, the meson masses Mb = 5.28GeV, Mk* = 0.89Ge\^, the decay constants 
fK* = 0.217GeV, /|, = O.lQGeV |32j, the central value of the CKM matrix elements 
\Vtd\ = 0.0075, \Vtb\ = 0.9992, = 0.9745, |Kfe| = 0.0033 and the B meson lifetime 

tbo = 1.536J9S {tb± = 1.671ps) [Hj. 

If we choose the CKM phase angle a{(f2) = 95° ^1], then our our numerical results are 
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given in TABLE. 1, where = Arg{A\\/ Aq) and 0_l = Arg{Aj^/ Aq). From the table we are 
convinced more with our power counting stated in chapter 2. We also find the polarization 
fraction R\\ ~ i?x and relative phase is around 2.6 for the former 3 channels. This is good 
news both for us and PQCD, since the current B — > (j)K* data |T71|^, which is also governed 
by the B — > K* form factors, suggest R\\ ~ R±, 0|| ~ 2.3 and 0_l ^ 2.5. These data are 
contrary from those rescattering effects jSHj and seem to support the evaluation of the relative 
strong phase in PQCD. 



TABLE. 1. Helicity amplitudes and relative phases 



Channel 


BR(10-^) 








(f)\\{rad) 


(f)±{rad) 




3.5 


0.78 


0.12 


0.10 


2.8 


2.8 


B+ K*+K*^ 


4.0 


0.75 


0.10 


0.15 


2.6 


2.4 


B- K*-K*° 


5.5 


0.88 


0.08 


0.04 


2.7 


3.0 


K*+K*~ 


0.22 


0.99 


0.005 


0.005 


4.1 


2.2 


BO ^ K*+K*- 


1.1 


0.99 


0.005 


0.003 


3.6 


1.9 



To test the contribution from different parts separately, we take B^{B^) x*^K*^ for 
example and classify the contributions into 4 kinds (see TABLE. 2.): (1) full contribution, 
(2) without annihilation nor nonfactorizable contributions, (3) without annihilation con- 
tributions, (4) without nonfactorizable contributions. Form the table we are convinced the 
annihilation contribution play an important role to the branching ratio. The annihilation di- 
agrams counteract with emission diagram severely, it even makes the branching ratio smaller 
when compared with the pure emission contribution. We also notice that contribution from 
the factorizable parts of the emission diagram is also bigger than the total branching ratio 
for B{B) — > As a result, if we change the form factor Aq from 0.4 to 0.32 as 

then our calculation give a branching ratio of 2.3 x 10^^. 
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TABLE. 2. Contribution from different parts: 
(1) full contribution, (2) without annihilation 
nor nonfactorizable contributions, (3) without 
annihilation contributions, (4) without nonfa 
-ctorizable contributions. 



Class 


BR(10-^) 








0|| (rac?) 


(f)x_{rad) 


(1) 


3.5 


0.78 


0.12 


0.10 


2.8 


2.8 


(2) 


4.4 


0.94 


0.03 


0.03 


TT 


TT 


(3) 


3.8 


0.86 


0.07 


0.07 


3.3 


3.3 


(4) 


4.1 


0.87 


0.07 


0.07 


2.5 


2.6 



For 5+ j^*+j^*o g^j^^ j^- _^ K*~K*^, it is similar to do so, we find annihilation 
diagrams contribute 7% and 31% to the total branching ratio respectively. If we take a 
smaller form factor and immediately get these values shown in TABLE. 3. We have to say 
these results are rather roughly because no precision wave function for K* have been given 
in that paper. 

TABLE. 3. The impact of a smaller form 
factor Aq = 0.32 on different channels 
Decay Channel BR(10-^) \Ao\^ l^yp |A^|2 

B^{B^) ^ K*^K*^ 2.3 0.67 0.18 0.15 

^ K*^K*^{K*^) 3.3 0.75 0.13 0.12 

To extract the CPV parameter of 5+ x*~^K*° and B^ K*^K*^, we can rewrite 
the helicity amplitude in ()18pi9|l as a function of the CKM phase angle a: 

Mjj = v:th-v;ph 

= V:Th{1 + ZHe'^'^^^"^) (32) 

Mh = VuTu-VtPH 

= KTi/(l + (33) 

where Zh = \V* /V*\\Ph /Th\, and 5 is the relative strong phase between tree(T) and 
penguin(P) diagrams. Here in PQCD approach, the strong phase comes from the non- 
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factorizable diagrams and annihilation diagrams. This can be seen from Eqs. ()75llOOI101|) . 
where the modified Bessel function has an imaginary part. This is different from FA ^ 
and Beneke-Buchalla-Neubert-Sachrajda(BBNS) |3j approaches. In that approaches, anni- 
hilation diagrams are not taken into account, strong phases mainly come from the so-called 
Bander-Silverman-Soni mechanism [33]. As shown in [7j, these effects are in fact next-to- 
leading-order(as suppressed) elements and can be neglect in PQCD approach. We give the 
averaged branching ratio of B"^ K*^K*^{K*^) as a function of a in Fig.4. 

5.00 

^ 4.95 I- 



CD 




4.60 - 



0.0 0.1 



0.8 0.9 

oJn 



Figure 4: Average branching ratios of K*'^K*'^{K*°) as a function of a. 



Using Eqs.(j22E2I), the direct CP violating parameter is 



A dir 
^CP 



|M+|'- 




|M+ 







—2sina (Tlsin6L + 2T^sin6N + ^T^sindr) 
Tl{l + ZI + 2ZLCOsacos6L) + 2 Ei=w,T (1 + + 2ZiCosacos6i) ' ^ ' 

Since the transverse polarization is twice of freedom when comparing with longitudinal one, 
the factor before Tn and Tt is twice as Tl. If we choose a as 95°, then the direct CP 
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asymmetry A'^'p for these channels are: 

A'^Sp{B\B^) K*^K*°) = 0, (35) 
A'^^''p{B^ ^ K*^K*%K*^)) = -15%, (36) 
A'^}.{B%B'^) K*+K*~) = -65%. (37) 

We notice the CP asymmetry of B^{B^) j{*oj(*o jg ^ero, since only pure penguin con- 
tribution in this channel. The CP asymmetry of B^ K*^K*^{K*^) is relatively small 
but large in B^{B^) K*~^K*~ , this is consistent with PQCD prediction. Using the power 
counting rules we stated in section. 2, the former channel is penguin dominated while the 
latter one is tree dominated, then from the definition of Zh we can easily deduce a big 
Zh for the former channel and a small Z^ for the latter one, so we can forecast a similar 
conclusion using Eqs. ()34j) without any calculation. We also notice the CP asymmetry for 
these channels are sensitive to a, hence we put Af^p as a function of a in Fig.5. 




.10 I ' I I I I I I I I I I I I I I I I I I I 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

a/n 



Figure 5: A-^^ as a function of a. (a) B°{B°) K*+K*-, (b) B± ^ K*^K*°. 

For the B^{B^) K*^K*~ decays, it is hard to distinguish B^ and -B", we can use the 
value given in TABLE. 1 to get an average branching ratio of 6.3 x 10^^. If we let CKM angle 
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a as a free parameter, then the evaluation of averaged branching ratio is shown in Fig. 6. 




1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 6: Average branching ratios of B^{B") K*+K* as a function of a. 

When it refered to CP asymmetry of B^{B^) decays, it is more comphcated due to the 
— B^ mixing. The CP asymmetry is time dependent: 

Acpit) ~ AcpCos{Amt) + a,+,i sin{/S.mt) (38) 

where Am is the mass difference of the two mass eigenstates of neutral B mesons. The 
direct CP violation parameter is already defined in Eg. ()34j) . while the mixing-related 
CP violation parameter is defined as 

—2Im{Xcp) 
= l + |AcP 

where 



2 



(39) 



In these two channels, over 90% of the branching ratios are composed of longitudinal fraction, 
so we can neglect the transverse contribution and use Eqs ()20l21|) to derive 
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People usually believe the penguin diagram contribution have been suppressed when com- 
paring with the tree contribution, i.e. in B{B) tt+tt" decay, Wilson coefficients of penguin 
diagram are loop suppressed, people believe Z <ti 1, \cp — exp[2ia], a^+e' = —sin2a, and 
~ 0, then it is easy to extract sin2a through the measurement of CPV. However, Z is 
not very small and a^+e' is not a simple function of —sin2a even in B{B) tt+tt^ [8], hence 
we couldn't get an exact a and this is called penguin pollution. In our channel, it is similar. 
After taking into account of the Wilson coefficient and the penguin enhancement we have 
stated, we get a Z^ as large as 0.80 and 5l = 2.35. We put a^+^z as a function of a in Fig. 7 
and we can see there isn't a simple relationship between fle+e' and —sin2a. 

1.2 I 1 




-1,2 I I I I I I I I I I I I I I I I I I I I 
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Figure 7: CP violation parameters a^+e' of B^{B°) —> K*+K*~ as a function of a 

If we integrate the time variable t of Eq. (jnH|) . we will get the total CP asymmetry as 

AcP = r^'^cP + T^«^+^' (42) 
1 + 1 + 

with X = Am/r ~ 0.771 for the B^ — B^ mixing in SM The integrated CP asymmetries 
for B^{B^) — > K*^K*~ are shown in Fig.8. As for B^{B°) K*^K*'^, there is only penguin 
contribution in this decay, direct CP is zero in Eqs.(jnH). The weak phase of penguin VtbV^*^ 
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is cancelled by the 5° — 5° mixing phase V^lVtd, so Xcp (see Egs- fHUj) ) is real here and 
a^+e' = 0. In fact, to the next leading order, there is a small up quark and charm quark 
penguin contribution which may give a small direct and mixing CP. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

aJn 



Figure 8: Acp of B{B) K*+K*- as a function of a 

When the PQCD formalism is extended to 0{al), the hard scales can be determined 
more precisely and the scale independence of our predictions will be improved. Before this 
calculation is carried out, we consider the hard scales t located between 0.75 — 1.25 times 
the invariant masses of the internal particles. For example, we take te{see Eqs. (j77j) ) in the 
following range, 

max(0.75Vi;MB, l/h, I/63) < < max{1.25^MB, l/h, I/63) , 
max(0.75v/x7MB, 1/61, 1/63) < t^e^ < max(1.25Vx7Mij, I/61, 1/63), (43) 

in order to estimate the O(a^) corrections. Then we can obtain the value area of the 
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branching ratio for the penguin dominated modes 

BR{B'^{B'^) = (3.5^^:3) X 10-^ (44) 

BR{B^ K*^K*°{K*°)) = {4.8tli) x 10-^ (45) 

which is sensitive to the change of t, so we can estimate that the next to leading order cor- 
rections will give about 20% contribution. The ratios Rq, R\\ and R± are not very sensitive 
to the variation of t, since the main contribution Flca, Fn^a and Ftea vary similarly when 
we conduct such changes on the maximum energy scale t. For the tree dominated modes 
B^{B^) —>■ K*^K*~ , the evaluation of the Wilson coefficients Ci and C2 is slow, hence it is 
not sensitive to the scale t, we changed the parameter uoh of the B meson wave function from 
0.32 to 0.48 and found the absolute value for each integration become larger when uo^ = 0.32 
and smaller when Uh = 0.48, then the value area for these decays are 

Si?(5°(fiO) ^ K*+K*-) = (6.41°;^) x 10"^ (46) 

If we compare our predictions with generalized FA j2] 

BR{B%B^) K*^K*°) = 3.5 x 10~^ (47) 
BR{B^ K*^K*^{K*°)) = 3.7 x 10"^ (48) 

with N^ = 3 and QCDF ^ 

BR{B^ K*^K*^) = 3.2 x 10"^ (49) 
BR{B- K*-K*^) = 3.4 x 10"^ (50) 

with the form factor from Light cone sum rules(LCSR) |22t I35j. we can easily find out their 
predictions for B^ K*^K*^{K*^) are a little smaller than our's, and neither of them 
gives the branching ratios and CPV parameters of B^{B^) K*^K*^, because annihilation 
diagrams can not be calculated in FA or QCDF in principle, while in PQCD all diagrams 
are calculated strictly. If we drop the annihilation contributions, we can get a similar results 
(3.7 X 10^'') with them and they are already shown in TABLE. 2. Current experiments jT3] 
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only give the upper limit for these decays 



^ K*^K*^) < 2.2 X 10-^ 
BR{B^ K*+K*^) < 7.1 X 10-^ 
BR{B K*+K*-) < 1.41 X 10"^ 



(53) 



(51) 



(52) 



and future experiments are expected. 

4 Summary 

In this paper we have predicted the branching ratios, polarization fraction and CP asym- 
metries of S — > K*K* modes using PQCD theorem in SM. We perform all leading diagrams, 
including both emission and annihilation diagrams, with up to twist-3 wave functions. The 
predicted branching ratios are compared with experiments and values from other approaches. 
We analyze the contribution from each parts for each decay channel, and found the anni- 
hilation diagrams is not very small to be neglected, then we present the dependence of the 
CP asymmetry and branching ratios on the CKM angle a. We also discussed the potential 
impact of a smaller form factor Aq in our paper. 
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A Factorization formulas 



The factorizable amplitudes are written as 

"1 



'Le4 



Te4 



ra3(4) 



dxidx2 



bidbib2db2^B{xi,bi) {[(1 + a;2)$if. (xa) 

+ rK*{l - 2x2){^'k-'{x3) + $^»(X2))] Ee4{t'^^^)he{xi,X2, 61,62) 
+ [2rK*$^*(x2) +r|*0K*(x2)] Ee4{tf^)he{x2,Xi,b2,bi)} , 



dxidx2 



bidbib2db2^B{xi, bi) 

Xrx* {[^K*{X2) +'2rK*^'^^{x2) +rK^X2{^K*{x2) - $^*(X2))] 

x^e4(tW)/ie(xi,a:;2,6i,62) 

+rK4^l,{x2) + ^'k'ix2)]Ee4{tf^)he{x2,Xi,b2,bi)} , 
rl poo 

IBttCfM^ / dxidx2 / bidbib2db2^B{xi,bi) 
Jo Jo 

XTK^ {[<^l.{x2)+2rK^<^l.{x2)+rK*X2{^''K'{x2) ' ^1*{X2))] 
XEe4t^^^)he{xi,X2,bi,b2) 

+rK4^l4x2)+^%.i^2)]Ee4itf)Kix2,Xi,b2,bi)} , 



327rCi;'M^rj^., / dx2dx2, f b2db2b^db^ 
Jo Jo 

X [(2 - X2){n4^2)n4^3) - n*{^2m*{^3)) 

+X2{<^UX2)<I>UXS) - <^UX2)<I>UXS))] 
X^ll(2)(^ii])^a(l-a;2, 1-^3,62, 63), 

327rCirM|r2, 



b2db2bzdbz 



dx2dx3 

X [(2 - X2){n4^2)n4^3) - $^(X2)$^(X3)) 

+X2{n4x2)n4xs) - ^Ux2)n4xs))] 

d) 

,3(4)1 



x^l?(4)(4Va(l - 2:2, 1 - a;3, 62, 63), 
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/•I POO 

= -327rCFM|r2. / c^xa^Xs / h2dh2hdh 
Jo Jo 

X [(2 - X2)in*i^2)n''i^3) - n*i^2m*i^3)) 

XEiS{C)ha{l -X2A- X3, 62, &3), (59) 
/•I /-oo 

Jo Jo 

X [(1 - X2)($|*(X2) + $*fe*(x2))$fc*(x3) + 2$fe.(x2)$^*(x3)] 

XEif{i^J)hail -X2,l- Xs, 62, &3), (60) 

^nIg = ^'^t^CfMqTk* / dx2dx3 / b2db2hdb3 

Jo Jo 

X{^1.{X2) i^Uxs) - ^U^s)) E^£{t'^])ha{l -X2,l- X3, 62, bs), (61) 

= 24^,, (62) 

p1 poo 

FtIs = ^'inCpMlrl, / dx2dx3 / b2db2bsdbs 

Jo Jo 

X [(1 - xs){rMrM - (t>Ux2WAx3)) 

+ (1 + X3)(0^(X2)0^,(X3) - 0^(0:2)^ (0:3))] 

x£;2(iW)/i„(x3,X2, 63,^2), (63) 

4^5 = -423- (64) 

The last expression of the factorizable amphtudes -4„5 doesn't really mean it equal to f!^^^ 
but with the evolution factor E^^^ replaced by eI^]^^ and plus a factor —1 in the beginning. 
Other amplitudes which you can not find in the upper formulas must be equal to zero. 

The factors E{t) contain the evolution from the W boson mass to the hard scales t in 
the Wilson coefficients a{t), and from t to the factorization scale 1/6 in the Sudakov factors 
S{t): 

E^Sit) = as{t)a^^\t)SB{t)SK*{t) , 

El^it) = as{t)af\t)SK.{t)SK,{t). (65) 
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The Wilson coefficients a in the above formulas are given by 



a? 




a?) 










( C-i 






"6 





+ leg ( Cg + 
+ ^eg ( Clo + 



3 

2 ^ 

3 

2 ^ 



^' 
No. 

Cs 

CV 
AT. 



(66) 
(67) 
(68) 
(69) 
(70) 
(71) 



kr resummation of large logarithmic corrections to the B, K* and K* meson distribution 
amplitudes lead to the exponentials Sb, Sk* and Sk*, respectively. 



Ssit) = exp 
SK*{t) — exp 
Sk* (t) — exp 



Ji/bi 

-s(x2P+, 62) -s{{l- x,)P^, 62) - 2 f ^7(«.(/^')) 

^1/63 ^ 



(72) 



The variables 61, 62, and 63, conjugate to the parton transverse momenta kir, k2T, and 
ksT, represent the transverse extents of the B, K* and K* mesons, respectively. The quark 
anomalous dimension 7 = — as/vr and the so-called Sudakov factor s{Q, b) is expressed as 

li/b 1^' 



s(Q,b) 



|(27.-l-log2)+C.log^^«^(^'^ 
3 jj,' 



,'67 10 2 7^1 f as{fi') 

+ S Uf + -po log — 

'9 3 27 ^ 3 2 



TT 



1 ^ 

log — 



(73) 



The above Sudakov exponentials decrease fast in the large b region, such that the B 
K*K* hard amplitudes remain sufficiently perturbative in the end-point region. 
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The hard functions /I's are 



he{Xi,X2MM) = Kq {y/x^MBh) St{x2) 

X [e{h - 62)^0 iV^MBh) Jo (v^Mfifoa) 

+e{b2 - bi)Ko (v^Mb62) /o iV^MBbi)] , (74) 

X e{b2-bs)H!,'\V^2MBb2)Jo{V^2MBbs) 

+d{b, - 62)//^'^ (v^Mb&3) ^0 (v^Mb62)] . (75) 

We have proposed the parametrization for the evolution function St{x) from threshold 
resummation[T^ ■ 

St{x) = ^'^rfi+t)'^ ^^^^ " ^^^l' • ^^^^ 

where the parameter c is chosen as c = 0.4 for the B K*K* decays. This factor modifies 
the end-point behavior of the meson distribution amplitudes, rendering them vanish faster 
at X — s> 0. Threshold resummation for nonfactorizable diagrams is weaker and negligible. 
Kq,Iq, Hq and Jq are the Bessel functions. 

The hard scales t are chosen as the maxima of the virtualities of the internal particles 
involved in the hard amplitudes, including l/bf. 

= max(v/i;MB, 1/61, 1/62) , 
tf) = max(V^MB, 1/61, 1/62) , 



ial(u) = max(v/r^MB, I/62, 1/63) , 
i« = max(yi^MB, 1/62, 1/63) • (77) 
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B Nonfactorization formulas 



The nonfactorizable amplitudes depending on kinematic variables of all the three mesons, 
are written as 

pi poo 

MLe3 = IQuCpMl^W, / d[x] / bidhih2dh2^B{XiM) 

Jo Jo 

xE'^,{t'i^)h'^^\x,,X2,xs,h,b2)}, (78) 

p1 poo 

MNe3 = I^'kCfMI^JWcTk* / d[x] / hidhih2dh2<^B{xiM) 

Jo Jo 

{X3$^.(X2)($^.(X3)-$^.(X3))£;:3(4'^)4'H^1,^2,X3,&1,&2) 

+ [(1 - a;3)$J(x2)($^*(x3) - $;(X3)) + 2rfc*(l + X2 - x^) 

x($^,(x2)*^.(a;3) -*^^.(x2)*^^.(x3))]£;:3(^f )4'n^i,^2,X3,6i,62)} , (79) 

MTe3 = 327rCFM|^2iVerx* / d[x] / hidhih2dh2^ b{xiM) 

Jo Jo 

{x3$^*(x2)($^.(x3)-$^*(x3))£;:3(tfVi'^(^l,^2,X3, 61,62) 

+ [(1 - a;3)$J(x2)($^(x3) - $;(a:3)) + 2ru*{l + X2 - x^) 

X($^.(X2)$^.(X3) -*^,(x2)*^.(x3))]£;:3(^f (^l,^2,X3,6i,62)} , (80) 

MLe5 = leTrC^Mlv^r^. / d[x] / hidh^h2dh2^B{XiM) 

Jo Jo 

{ [X3^K*{X2){^K*{X3) - ^4,{X3)) + rK*{x2 + X3) ($^, (^2)*^* (^^3) + ^K*ix2)^K'ix3)) 
+rK*(x2 - Xs)(^U^2)^K*(xs) + ¥j,4x2)^U^3))] KSf)hf (^Xl,X2,X3,0i,02j 
+ [(-1 + X,)^k*{x2Wk*{^3) + ^1{X,)) 

+rx*(-l - X2 +X3)($^*(X2)$^*(X3) - ¥j^.{x2)^K*{x3)) 

+rK*(-l + X2 + X3)($^,(a;2)$*^.(a;3) -$*K*(a;2)$^.(a;3))] 

XE'^,{tf)hf{x,,X2,X3MM)} , (81) 
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/•l poo 

MNe5 = IdTrCpMly^^rk* / d[x] / hdhb2db2^Bixi,bi) 

Jo Jo 

{ [X2($^(X2) - $^.(X2))$J(X3) + rfc.(x2 + Xa)*^- (2^2)$!. (2^3)] 
xKb{tf)hf{xi,X2, X3,bi, 62) 

+ [X2($^(X2) - ^*{x2))^l.ix3) + rif*(l + X2 - X3)$[«(x2)$fe«(x3)] 

xE'^,{t^J'^)h^P{x^,X2,Xs,b^,b2)} , (82) 



MTe5 = 327rCFM|y27Verfc. / o?[a;] / b^db^b2db2^B{xM 

Jo Jo 

{ [X2{n*{x2) - <f^.(a;2))<l'J(x3) + r,.(a;2 - X3)^l{x2)^l{x3)\ 
xK5i't'd^)hd\xi,X2, X3,bi, 62) 

+ [x2($fe.(x2) - $^*(x2))$fc.(x3) + rif*(-l + a:2 + X3)$fe*(x2)$fe*(a;3)] 

X^e5(4'^)4'^(a^l,^2,:r3,&l,&2)} , (83) 

/I /-co 
/ b,db,b2db2^B{Xl,bi) 

{[{-1 + X3)<^>k*{x2)<^>k'{x3) 

+4* {-^H'{X2)H4^3) + {X2 + X3)m,{x2)^,{x3) + ^1,{X2)^1,{X3)) 

+ {X2 - X3){n*ix2)^Uxs) + ^UX2WA^S)))] 

y~E%^{tf)hf,{x,,X2.X3MM) 
+ [(1 -X2)$fe*(x2)$fe*(x3) 

+rl ((2 -X2- X3){<^UX2)<^UX3) + ^\.{X2)^\.{X3)) 
+ {X2 - X3){^l*{x2)^i,{x3) + ¥^,{X2)^1,{X3)))] 

xE%{t^P)hfJ{xr,X2,X3A,b2)} , (84) 

/» 1 /'OO 

•^!^i3(4) = 167rCfM|y2iVe4. / / bidb,b2db2^B{Xl,b,) 

{ [2{^'k.(x2)^'k*(x3) - ^'k*(x2)^K^(x3)) + (x2 + 2:3)$^. (^2)$^. (^s)] 

x^2(4) (^1, 2:2, a;3, bi, 62) 

+ (2 - X2 - X3)$^* (X2)$^* ix3)E%) iif)hfd {xi, X2, X3, 61, 62) } , (85) 
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pi poo 

^Ta3{4) = ^2T:CpMl.j2N,rl. d[x] hdhh^dh^eixM 

{ [2($^.(X2)$^.(X3) - $^.(X2)$^.(X3)) + (X2 - Xg)*^. (Xs)*^. (Xg)] 
X ^if {tf ) hfd {xi,X2,X3,h,b2) 

+ {Xs - X2)^l. {X2)^l. {^S)E%) {tf)h% (Xi, X2,Xs,h, 62) } , (86) 
pi poo 

■^il = 167rC^M|V27VerK* / d[x] / hdhh2dh2^ b{x^M) 

Jo Jo 

{[(1 +X2)($*;^.(X2) - ^%,{X2))^K*{X3) + {1+Xs)^k4x2){^%*{x3) - ^K'i^s))] 

XE^^'{tf)hf,{x,,X2,X3MM) 

+ [(1 - X2){^K*{X2) - ^K'{^2))^K*{xz) + (1 - Xs)*^* (2^2) (*^« (2:3) - ^K^ix^))] 

xE^^'{tf)hf^(x,,X2,xsM,b2)} , (87) 

pi poo 

M%i^ = 167^C^M|^/27VerK* / o?[x] / h^dh^h2dh2^ b{x^M) 

Jo Jo 

{[(1 +X2)($^*(X2) - ^1,{x2Wk,{x3) + (1 +X3)$^*(X2)($^.(X3) - ^^.(Xa))] 

X ^if ' (4^^ ) ^fd {xi,X2,X:iMM) 

+ [(1 - X2)($^.(X2) - $^.(X2))$^.(X3) + (1 - X3)$^,(X2)($^.(X3) " ^1c*{x,))\ 

xE^f/{tf)hf^{x^,X2,xsM,b2)} , (88) 
M?^, = 2M^il. (89) 

/»! poo 

M^l = WttCfM^sV^c d[x] hdhb2db2^B{xi,h) 

Jo Jo 

{[{-I + X2)^k*{x2)^k*{x3) 

+rl' {-^^1,{X2M,{X3) + {X2 + X3)m,{x2)^,{x3) + {^2)^1' i^^)) 
+ (x3-X2)m.(x2)^Uxs)+^l(x2)n.(xs)))] 

X Eif {tf^ ) hf^ {xi,X2,x^MM) 

+ [{1 - X^)^k'{x2)^k'{xz) 

+rl, ((2 -X2- X3)($^(X2)$^(X3) + $*fc.(x2)$*fc.(x3)) 

+ (X3 - X2)(<f^,(a;2)$*.(a:3) + ¥^.{x2WA^^)))] 

xE^£'{tf)hf^{x,,X2,x^MM)] , (90) 

A^!vl = M'^il (91) 

A^?^ = --MSs (92) 



27 



p1 poo 

Ml = -levrC^M^v^ / d[x] hdhb2db2^B{xi,bi) 

Jo Jo 

{ [X2'^k* ix2)<^k* ix3) + rl (2($|. ix2)<^l* (Xg) - (x2)$fc. (xg)) 

+ (X2 + X3){^l,{x2)^l,ixs) + $*fe,(a;2)$t.(x3)) 

+ (X2 - X3)m.{x2)^Uxs) + $|.(X2)$^(X3)))] 

Xi?(f(4^))4^i(Xi,X2,X3,6i,62) 

+ [-X3<^>k* {X2)<^k* {Xs) + {-{X2 + X3) {^l. {x2Wk* (^3) + (X2)$fe« (X3)) 
+ (X2 - X3)($^*(x2)$*fe*(x3) + $|.(X2)$^*(X3)))] 

xE^I'{tf)hf{xuX2,X3MM)} , (93) 

/I /-oo 
d[x\ / bidbib2db2^B{xi,bi) 

{ [2($^,(a;2)$^.(x3) - $^.(0:2)$^ (0:3)) + (2 - X2 - 0:3)$^. (0:2)$^. (0:3)] 
x£;S'(4^))/ig(a;i,X2,X3,6i,62) 

+ (X2 +X3)$^.(x2)$^.(x3)i?lJ(l)(f V^'i(xi,X2,X3,6l,62)} , (94) 

/I poo 
d[x] / bidbib2db2^B{xi,bi) 

{[2($^,(a;2)«f^(x3) - $^ (0:2)$^^. (0:3)) + (0:2 -0:3)$^. (0:2)$^. (0:3)] 
xS(f(4^))4^,)(a:i, 0:2,0:3,61,62) 

+ {xs-X2)^l4^2)^l4^s)E%fif)hf^{x,,X2,X3,b,,b2)} , (95) 

A^(.*^ _ -167rCj.M|V27Ve / d[x] / 61^6162(^62*5(2:1,61) 

{[a:3$fc*(a:2)$fc*(x3) +r^* (2($^.(x2)$^*(a:3) - $*fe*(a:2)$*fe*(x3)) 

+ {X2 + Xs){<!>l4x2)n4x3) + <^UX2WAXS)) 
+ {X3 - X2)m.{x2)^l{x3) + ¥^4X2)^4^3)))] 

XE^^'{tf)h^l{x,,X2,X3MM) 

+ [-X2^k'{x2)^k'{x^) + 4* {-{x2+X3){¥k*{x2M*{x^) + ^\*{x2Wk*{x^)) 
+ {X3 - a:2)($^(x2)$*fe*(x3) + $i.(x2)$^(x3)))] 

X^i?'(tf )4'i(^l>^2,X3,6i,62)} , (96) 

= -Mi^i, (97) 
= -MS- (98) 
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The expressions of the nonfactorizable amphtudes Aina and AineA are the same as Ai^Ha3 ^^"^ 
A4^g3 but with the evolution factors E^^' and E^f' replaced by E^^' and E^f', respectively. 
The evolution factors are given by 



E^'{t) = a,{t)ar{t)S{t')\,,=,,, 



{<!)'( 



(99) 



with the Sudakov factor 5" = SbSk*Sk*- The Wilson coefficients a appearing in the above 
formulas are 



(I , — 



9i 



The hard functions j — ^ and 2, with d stand for emission and / stand for annihi- 
lation, are written as 



h^fg{xi,X2,X3,bi,bs) 



[0{bi - bs)Ko{DMBb,)IoiDMBbs) 

+9ib3 bi)KoiDMBb3)Io{DMBbi 
KoiD^MBbs), 



2 -"O ' 



D] > 0, 
D^j < 0. 



\D]\MBb^), 
e{b, - bs)Hl,'^ {FgMBb,) Jo (F.MBbs) 
+e{bs - (F.MBbs) Jo (F.MBbs) 



X 



KoiFjgMBbi), ; 



Fj.lMBbi 



(100) 



(101) 
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with the variables, 

= X1X2 , 

Dl = X2{Xi -X3) , 

Dl = -^2(1 -xi- X3) , 

Fd{u) = (1 - 2:2) (1 - X3) , 

Fld{u) = 1-X2(x3-Xi) , 

= X2a;3 , 

i^l' = 1-(1-X2)(l-Xi-X3) , 

Fl = X2 (xi - X3) . (102) 
The hard scales t^^^ are chosen as 

= max (^DMb, .^/^Ms, 1/61, 1/63^ , 

if = max (^L>Mb, y^Mb, 1/61, 1/63) , 
4;)' = max (f,Mb,.^/^Mb, 1/61, 1/63) , 

4f = max (f,Mb, .^/^Mb, 1/61, 1/63) ■ (103) 
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